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ABSTRACT 


An analytical model was constructed to predict the 
response Of tne ball obturator in a ball obturated, spin 
stabilized tubular projectile (BOT). Toxques acting on 
the obturator were considered to be due to fluid and sliding 
Friction forces. An experimental apparatus was designed, 
and built to simulate the BOT. Data obtained were the times 
required for the obturator to nutate through a known angle 
Over a range of projectile angular velocities. Utilizing 
published coefficient of sliding friction information and a 
linear approximation of the exact solution, the trend in the 
data to go through a minimum value was duplicated. A least 
squares fit of the linear approximation to the data allowed 
fiewrormulation of a new sliding friction coefficient, u.. 
mermmioetiis value Of t., a plot of the least squares fit, 
values from the exact solution and the experimental data 
was constructed. The correlation is extremely good so that 
confidence in the predictive capabilities of the model was 


gained. 
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NOMENCLATURE 


(Only primary symbols are listed. Intermediate quantities 


are defined in the text.) 


A major maSS moment in inertia, N-m-s 

ie minor maSS moment of inertia, N-m-s 

G applied torque, N-m 

h width of gap between ball and projectile, m 
M magnitude of applied torque, N-m 

ro radius of hole through projectile, m 

ie radius of hole through ball, m 

Ro radius of ball, m 

8 SHR Na De/ Re) 


8,o,W Euler angles (see Fig. 7) 


A C/E eal 

Us fluid viscosity (Pa-s) 

He coefficient of sliding friction, dimensionless 
u) angular velocity 

Subscrivots 

er refers to condition when 0=0.4r=B8s 

O initial value 

D projectile 

i refers to relative angular velocity 

S Sphere (ball obturator) 


eb 





cay, Z refers to inertial coordinate system 


ime, 3 refers to body-fixed coordinate system 
Superscripts 
. denotes differential with respect to time 
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L. INTRODUCTION 


A. BACKGROUND 

The tubular projectile is not a new concept. As early as 
1858, the usefulness of tubular projectiles was recognized for 
Specific purposes. In that year, Joseph Whitworth (Whitworth 
Threads) wrote about and included an illustration of the tubular 
projectile in the section on Rifled Firearms of his Miscellaneous 
Papers on Mechanical Subjects [1]*. The projectile pictured 
was hexagonal with a circular hole. Whitworth noted its par- 


ticular "...effectiveness in perforating elastic materials 
which prevented them from closing up." The projectile utilized 
a wooden sabot. Whitworth also said that the tubular projectile 
penetrated deeper into masonry than any he was acquainted with. 
The Krnka-Hebler projectile was the next projectile refer- 
enced chronologically [2]. According to Charters and Thomas, 
this projectile was reported in the Allgemeine Schweitzerische 
Militarzeitung as having been so successful that the Ordnance 
Department in the United States carried out firing tests of 
caliber .30 tubular projectiles in 1894. A description of 
those test firings and an updating of the results was performed 
by I.E. Segal and included in Ref. 2. In the 1894 tests, a 
vulcanized fiber sabot was used to push the projectile out the 


parrel. 


1 : : 
Numbers appearing in brackets refer to the list of references. 
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The purpose of the 1894 test was to compare trajectories 
of standard and tubular projectiles. This was done by com- 
paring the vertical drop on target at a given range. 

Segal's report in Ref. 2 indicates that the drag co- 
efficients computed from the 1894 results agree closely with 
those of Charters and Thomas [2], even though their report 
was published fifty years later (1944) and their experimental 


apparatus wasS more sophisticated. 


Be RECENT DEVELOPMENTAL STUDIES 

The results obtained by independent researchers as to the 
value of the tubular projectile as an alternative to the con- 
ventional round are contradictory. It seems as if there was 
as much research devoted to discrediting the tubular pro- 
jectile as there was to investigating itS merits. Frank and 
McLaughlin [3] have accumulated a great deal of data from 
Various sources in an attempt to "objectively compare" the 
merits of the tubular and "conventional shapes". The authors 
concluded that the tubular has no particular advantages over 
"well designed" conventional projectile shapes. 

Their findings disagree with most of those discussed in 
the following sections. 

me Range Tests 

Winchenbach, Daniel and Edgar [4] conducted range 

tests of six configurations of tubular projectiles and con- 


cluded that the drag coefficients were significantly lower 
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than the standard High Explosive Incendiary (HEI) projectile 
of the same caliber. Only projectiles of the same bore size 
were compared. The experimental models were constructed from 
standard 20-mm ammunition by boring various hole sizes to 
@emotin ". . . area ratios swe7A7eNe Ouslig,. OacwmeOeomand 1.30..." 
Where At is the throat area and A; is the inlet area of the 
projectile. The drag coefficients of the projectiles with 
area ratios of 0.8, 0.9 and 1.0 were less than half that of 
the projectile with area ratio of 0.7. The higher drag co- 
efficient was attributed to the underexpanded flow at the 
exit resulting in high base drag. 

Range tests conducted for the concept evaluation of 
the 20-mm tubular projectile for the Vulcan Weapons System [5] 
showed that the tubular projectile performed better than the 
Standard HEI round,with lower drag, shorter time of flight for 
a given range (30 percent shorter at 1000 meters, 40 percent 
shorter at 2000 meters) and similar dispersion characteristics 
Gaecarget. 

Recent tests at the Naval Weapons Center, China Lake, 
California of the Ball Obturated Tubular Projectile (BOT) 
have shown a definite advantage in that the low drag of the 
tubular projectile means slower retardation of velocity and 
hence a higher terminal kinetic energy [6]. (The BOT is the 
object of this investigation and will be described in detail 


in later sections of this thesis.) 


Ey 





a aelercetyinpact Tests 
The higher the energy delivered to the target, the 


greater the damage. Since kinetic energy is directly propor- 
tional to the first power of mass and to the square of the 
velocity, the reduced mass of the tubular projectile is com- 
pensated by a higher muzzle velocity and therefore a potential 
for higher kinetic energy. However, due to the improved drag 
characteristics, the reduction in mass to obtain the higher 
muzzle velocities may not be necessary, depending upon the 
imposed performance criteria. 

Target impact studies have been carried out with 
various projectile configurations and target types. 
Rhethorst, et.al. [7] conducted impact studies of 7.62-mm 
tubular projectiles on helmets. Tests showed that even with 
the same energy, the tubular projectile penetrated further. 

Kitchen and Keeser [8] conducted studies for the 
Air Force on the impact effectiveness of tubular projectiles 
On Simulated aircraft fuel cells. These tests were conducted 
with steel and depleted uranium (DU) tubular projectiles and 
Standard 20-mm HEI projectiles. The projectiles were fired 
at double panels at varying degrees of obliquity. Of the 
forty-seven tests of the steel tubular projectiles which im- 
pacted the target, twenty breached the rear panel up to 
angles of 70°. The DU projectiles breached the rear panel 
even at angles of g5° and fires were started in three of the 


five DU tests. The standard M56 HEI projectiles failed to 
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breach the rear panel in each of the five firings even though 
severe damage to the front panel and two fires occurred. 

Brunsvold and Kalivretenos [9] conducted a program 
to test the effectiveness of the tubular projectile against 
a Simulated cruise missile warhead. Tests were conducted 
with 20 and 40-mm rounds. Results of those tests are 
classified, and therefore not presented here. It was noted 
[9] that due to the improved stability characteristics of 
the tubular projectile, only half the spin rate need be im- 
parted by the rifling in the barrel to obtain the same sta- 
bility as a standard projectile, which indicates less barrel 
wear. 

3. Weapons System Capability 

For the tubular projectile to be practical, it must 
be able to interface with existing weapons syStems. Refer- 
ence 5 1s a report of the 1978 tests conducted by the Army 
to evaluate the tubular projectile in the Vulcan weapons 
system. 

Because the nose of the tubular projectile is flat 
compared to the standard projectile (See Fig. 1), it"... 
did not lend itself to chambering in the weapon." Personnel 
safety required that the weapon be remotely operated and 
therefore only the surface-to-surface performance comparison 
was made. 

The dispersion on target improved for the tubular 


projectile from 700 to 2000 meters. The dispersion varied 


a9 











Figure l. 








ties 





20-mm Ball Obturated Tubular Projectile (BOT), 
components, and conventional 20-mm projectile. 
(Photograph courtesy of NWC, China Lake). 





Figure 2a. 


Standard 20-mm projectile in flight with 
detached bow shock. (Photograph courtesy 
of NWC, China Lake). 
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Figure 


Figure 





2b 


BOT with passage blocked resulting in 
detached bow shock. (Photograph 
courtesy of NWC, China Lake). 





ZC. 


BOT with no blockage in passage; only attached 
oblique shocks emanating from the lip of the 
projectile. (Photograph courtesy of NWC, China 
Lake). 
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on the average by only 0.2 mil over the entire range spectrum 
between the tubular and conventional HEI projectile. However, 
at 2000 meters, the tubular was better by 0.4 mil. The time 
Of flight to target was another impressive finding in Ref. S. 
The tubular projectile had a 30 percent shorter flight time 
at 1000 meters and a 40 percent shorter flight time at 2000 
meters than the conventional HEI projectile. 

Results of firings of 25-mm and 30-mm tubular pro- 
jectiles from the Oerlikin KBA and GAU-8/A respectively were 
reported in part in Ref. 3. Figures comparing Line-of-Sight 
Penetration Capability vs. Range for a conventional spinner 
configuration based on AR-2 shape, a finner configuration 
based on F1lO shape, and a tubular projectile were shown. The 
tubular projectile was shown to be inferior in both instances. 
The report does not indicate the degree of compatibility the 
tubular configuration has with the guns used. 

4. Review Of Internal Shock Wave Considerations 

From photographs of tubular projectiles in flight [4] 
and during wind tunnel testing [7], various shock patterns 
have been observed at the inlet region and in the wake. The 
photographs of Figs. 2a-c show the BOT and the conventional 
projectiles in flight. The sabots used in Refs. 4 and 7 were 
of the puSher type similar to those illustrated in Fig. 3. 
The detached bow shock seen in the photograph of the standard 
Beenvectile (Fig. 2a) 1s Similar to that of the tubular pro- 


jectile shown in Fig. 2b with passage closed (and in Ref. 4 


pe 
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With the sabot attached). With no blockage of the internal 
passage (Fig. 2c) there was no detached bow shock, with only 
attached oblique shocks emanating from the lip of the pro- 
jectile. This latter case allows the possibility of four 
other internal flow/shock configurations. 
(1) A normal shock standing at the entrance to the pro- 
jectile; 
(2) A normal shock standing at some intermediate position 
in the channel; 
(3) A system of oblique shocks present in the channel; 
(4) The channel is devoid of shocks and the flow is shock- 
free through the channel. 
In each of the cases described, specific conditions 
must exist: 
a. The Detached Bow Shock 
When the projectile passage is closed, such as 
in the case of the attached sabot, the projectile acts as a 
blunt-nosed body with the resultant detached normal shock [10] 
(see Figs. 2a and b). From oblique shock theory, there is a 
maximum turning angle which the flow (Supersonic) can nego- 
tiate through an attached oblique shock. When this turning 
angle is exceeded, the shock becomes normal and detaches. 
There have been many studies as to the strength and location 
of this detached bow shock, but here the internal flow is of 
more interest and, therefore, those references are not in- 


Gidea in this discussion. 
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b. Normal Shock At The Entrance 
From one-dimensional frictionless flow theory, 
for a given initial Mach number, Mos the pressure ratio 
f@enosS a normal shock is given by: 
P (1 + yM?) 
2 1 


= 54s (1) 
~s (ac YM) 


for a perfect gas with constant specific heats. P is the 
pressure downstream of the normal shock, e the pressure up- 
stream of the shock, af the downstream Mach number and y the 
ratio of specific heats, Co/ Cy: This pressure ratio must 
exist in order for a normal shock to be formed. Kantrowitz 
and Donaldson [11] conducted experiments in which they con- 
cluded that the contraction-ratio in the entrance of a Super- 
sonic diffuser was important in reducing the kinetic energy 
losses. The studies done were on Supersonic diffusers with 
varying areas. The formation of a normal shock at the en- 
trance is determined by the throat area and initial Mach 
number such that the Mach number at the throat of the duct 

is exactly 1.0. The maximum value of the ratio of the capture 


(inlet) area, A to throat area, A,, required to allow the 


c?/ 


bow shock to be swallowed iS given in Ref. 9 as: 


v+1 
2(y-1 
me = 1 —— (1 + Y= My] ae (2) 
At M y+1 2 2 
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where M, may be obtained from tabulated values for normal 
shock waves for a given free stream Mach number, M or may 
1 


be calculated from 


a 


eel ine 
te = y= 1) ] (3) 


ZY og 
(yale) i : 


c. Normal Shock At An Intermediate Position 

The formation and position of a normal shock at 
an intermediate position in a constant area duct has been 
theorized in the presence and absence of viscous effects. 
In the absence of viscous effects, the normal shock tends 
to grow from the interaction of oblique shocks and will be 
discussed under "d" below. 

inmewe presence Of Exiction, the formation of 
the boundary layer along the duct walls reduces the effective 
peeemor the duct (12, 13, 14]. This reduction has the effect 
of decelerating a supersonic flow as in a diffuser. Ina 
constant area duct with friction, the location of a normal 
Shock in the duct may be determined by Fanno Line Flow Theory. 
The length to diameter ratio, friction factor and initial 
Mach number determine the range of back pressures for which 
the duct will be "choked", or a normal shock set up [10]. 

For a given Supersonic entrance Mach number, there 
exists a maximum value of fL/D; £, being the Reynolds number- 
dependent friction factor, L, the length of the duct and D, 


the diameter of the duct. This maximum value, is the value 
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fone ewhich the duct will remain unchoked. If the actual value 
of fL/D is less than the maximum, the flow will remain un- 
choked if the back pressure is sufficiently low at the duct 
exit. 

In the case of the tubular projectile, assuming 
a constant back pressure, L/D of (typically) 8.0 and Re in 
the range from approximately 4x10° to 6x10°, from Fig. 20.1 
of Ref. 15, f is between 0.013 and 0.0125. Using an average 
meme of f of 0.01275, 

£L/D = O7102 
From Ref. 10, Appendix E, this corresponds to M#1.4l. 
Iterating to obtain better values of f from a revised average 
Re would increase the Mach number at which the projectile 
would choke. After one iteration the Mach number increases 
to M= 1.44. As L/D decreases, Re increases and f decreases 
resulting in a lower value of the Mach number at which the 
flow chokes. 
d. System Of Oblique Shocks In The Duct 

The formation of a system of oblique shocks in a 
Guct or channel is the most difficult to explain and analyze. 
Ferri [{16, 17] has developed an analytical solution for the 
Shock structure utilizing the Method of Characteristics. 
The analysis becomes difficult if not impossible at the 
centerline of the duct because, according to Ferri, the oblique 
shock curves toward the centerline and becomes normal. The 


extent of the normal shock is dependent upon geometry and flow 
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conditions. The larger the wedge angle for a given Mach 
number, the stronger and larger the normal shock at the 
centerline. The constant area duct implies a zero internal 
deflection angle. In this instance, Ferri's analysis, which 
ignores viscous effects, implies that no normal shock is 
permed (Fig. 13 of Ref. [16]). With no wedge angle, there 
is no oblique shock formed. 

Brunsvold and Kalivretenos [9] have modified a 
Method of Characteristics computer code [18] developed for 
Supersonic inlets to account for the absence of a center-body 
on the tubular projectile. This modified code was then 
coupled with a subroutine to calculate the internal wave drag 
coefficient based on A, and used to obtain an optimum internal 
geometry for the tubular projectile which would give minimum 
internal wave drag. The computer code assumes inviscid, 
supersonic flow throughout. Reference 9 states that some 
inaccuracies are allowed in the computation due to the re- 
laxation of the convergence of the characteristics equations 
at the centerline, but that the inaccuracies do not have time 
to propagate to the walls of the projectile before the exit 
plane. Therefore, the inaccuracies do not affect the wall 
pressure distribution. 

Rhethorst, et.al. [19] has postulated another 
method for determining the shock structure and the resultant 
internal wave drag for a tubular projectile, using Prandtl- 


Meyer wave theory. But again the initial oblique shock 
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formation is dependent upon the wedge angle at the entrance 
and therefore does not apply to the constant area duct case 
where no internal wedge angle exists. Ferri [16, 17], 
Brunsvold and Kalivretenos [9] and Rhethorst, et.al. [19] 
have neglected boundary layer/shock wave interaction. 

Fejer, et.al., [14] and Waltrup and Billig [15] 
have conducted experiments utilizing constant area ducts. 
Fejer, et.al., performed experiments on constant area rec~ 
tangular ducts varying in lengths, whereas Waltrup and Billig 
conducted their experiments on cylindrical ducts of varying 
lengths and diameters. Other than the obvious geometry dif- 
ferences, the procedures were very similar. The advantage 
of the rectangular duct was the ability to visualize the shock 
patterns through a viewing port, whereas the curvature of the 
walls of the cylindrical duct precluded accurate visualization. 

In both experiments, oblique shock systems were 
formed in the ducts. The location of the shock system in the 
duct and the actual configuration of the shock system was con- 
trolled by the pressure ratio. The pressure ratios were con- 
trolled in different ways, but the effects were the same. As 
back pressure increased, a shock system formed at the exit of 
the duct. With further increases in back pressure the system 
of shocks moved upstream and out of the duct and a detached 
shock appeared. Shapiro [12] indicates that internal 


oblique shocks may be the remnants of normal shocks "with 
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bifurcated ends as a result of boundary layer separation." 
He states that the normal portion grows shorter as the 
boundary layer gets thicker until the normal portion dis- 
appears completely leaving only the bifurcated portion. He 
illustrates this in Fig. 28.27(c) of Ref. 12. This explains 
the movement of the oblique system toward the entrance of the 
duct in the experiments of Refs. 14 and 15. As back pressure 
increases, the flow downstream of an internal shock has more 
of a tendency to separate. Fejer also states that the 
presence of a turbulent boundary layer had a stabilizing 
effect on the shock system in the constant area rectangular 
duct. In the experiments involving the cylindrical duct, tne 
shock structure again moved from exit toward entrance with 
back pressure increases. They also report a breakdown in flow 
downstream of the shock system as indicated by relatively 
constant pressure movements in this region. 

e. No Shocks Present In The Duct 

It can be seen from previous discussion that 

given an initial Mach number and inlet stagnation pressure, 
a sufficiently high back pressure will cause a normal shock 
or system of oblique shocks to be formed. If the cross- 
sectional area of the duct is reduced in the direction of 
flow, aS in a Supersonic diffuser, the incident flow may 
cause an oblique compression wave to be formed. In the 


absence of the requisite back pressure or area reduction, 
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no explanations for internal shock formation have been found 
other than those due to viscous effects noted in "d" above. 
ff. Summary 

iM a COnStant area duct with supersonic inlet 
flow, there is a pressure rise caused by a reduction in the 
effective flow area by boundary layer formation due to 
friction. Depending upon the length-to-diameter ratio (L/D) 
Of the duct and the flow parameters, the increase in pressure 
ratio (or reduction in Mach number) required to cause forma- 
tion of a system of oblique shocks is less than that required 
to set up a normal shock at the duct entrance. This difference 
is explained by the boundary layer formation and/or interaction 
with the shock waves. 

In the case of the tubular projectile, the back 
pressure at the exit in flight is affected by the complex wake 
flow. There is no pressure difference in the surrounding 
medium at large distances from the projectile. The formation 
of shock waves inside the projectile with a constant cross 
section is due to viscous effects and depends upon Mach number 
and the length-to-diameter ratio. 

oe Gaboratory/Wind Tunnel Testing 
From all indications, the Weapons System Concept Team 
(WSCT) of ARRADCOM has done extensive wind tunnel 
testing on various tubular projectile configurations. The 
work of A. Flateau of the Weapons System Concept Office is 


referenced often with regard to results obtained, however no 
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publications from that office have been found. (The work 
referenced in other publications is usually of the "private 
communication" type.) Appendix D of Ref. 8 reports results 
Of a test conducted by the WSCT at the ARRADCOM facility. 
The facility utilizes a variable Mach number tunnel which 
may be used to simulate the deceleration of the projectile 
ier light. 

The results of Ref. 8 agree closely to those of 
Ref. 7, in that the drag coefficient rises almost instan- 
taneously at the Mach number at which the flow becomes 
choked. 

The Mach number at which this choking occurs is de- 
pendent upon geometry. For the case of inviscid flow, 
Brunsvold and Kalivretenos [9] give the following relation 
for the area ratio corresponding to the free stream Mach 
number at which the shock is "regurgitated" and the flow be- 


comes choked: 


c aL 2 vid 
= Se = =e See Se ene 
Ay M y+1 (2 Z ) (4) 


SeeeeeeRPOSE OF STUDY 

The results of tests and experiments previously mentioned 
show for the most part definite advantages to the tubular 
projectile over the conventional projectile primarily in the 


area of lower drag (approximately 1/2 to 1/3 of conventional [8]) 
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which allows a flatter trajectory and shorter time of flight 
for a given range [5]. With the emphasis on Anti-Missile 
Defense both from a Surface-to-Air and Air-to-Air standpoint, 
jeese features are definitely worth further study. 

The problem seen from the air platform is that of the 
discarding sabot or pusher. The solid [4, 5, 7] or split [8] 
disk is effective and reliable, however the possibility of 
ingestion in the aircraft engines prohibits their use. 
Rhethorst et.al. in Ref. 7 illustrates a number of sabot/ 
obturator designs for use with the tubular projectile, but 
these are also of the discarding type. They also investigated 
the possibility of using a consumable sabot [19] which would 
be burned up as the projectile was transiting the barrel and 
would be completely consumed as it exited the barrel. The 
most recent development has been the BOT which was designed 
memieeenwChina Lake, California [6]. The ball obturator has 
been bored with a hole the same diameter as the hole through 
the projectile. When loaded, the ball is supported inside 
the projectile such that the axis through the hole in the ball 
makes an angle of approximately 90° with the axis of the hole 
through the projectile (see Fig. 4). 

Gas pressure from the burning propellant holds the ball 
obturator fixed with respect to the projectile as the rifling 
in the barrel imparts spin to the projectile (as well as the 


ball). Upon exiting the barrel, the propellant gas pressure 
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1S released. At that time a complex force distribution acts 
en the ball to cause it to change position inside the pro- 
jectile in such a way as to align the holes. [Inertial forces 
then dominate to maintain this alignment so that the pro- 
jectile remains fully tubular. There are no separating parts 
and the opening process iS automatic. 

The purpose of this research is to analytically predict 
Mieemotion of the ball obturator as a function of time. in 
Support of this purpose, an experimental apparatus was de- 
Signed and built to simulate the motion of the obturator in 


Order to gain phySical evidence of the accuracy of the model. 
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It. ANALYSIS 


A. DEVELOPMENT AND SOLUTION OF THE EQUATIONS OF MOTION 

The ball obturator may be thought of as a rigid body of 
rotation with a system of coordinate axes fixed to the ball 
having its origin at the mass center of the ball (Fig. 5). 
These axes will be designated x, y, and z, where z is the 
axis through the hole in the ball. 

The motion of the ball must be described relative to the 
projectile, and for this purpose an inertial frame of refer- 
ence is defined with axes designated X, Y and Z with its 
Origin also at the mass center of the ball. For this analysis 
1t 1S assumed that the projectile has only a rotation about 
the @-axis and no translation or other rotation. Because only 
the motion of the ball with respect to the projectile is de- 
Sired at this point, the relative velocities are the quantities 
desired. The ball has only rotational components and thus the 
values of angular velocity with respect to the X and Y axes 
are also relative to the projectile, whereas the value with 
respect to the Z-axis must account for the projectile spin. 

By doing this, the motion may be described by Euler's 
Modified Equations of Motion for a rigid body of revolution 


about a fixed point. These are given in Thomson [20] as: 
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Aw + (C-A)w w =G (5a) 
a3 1 
A® + (A-C)w = G. (5b) 
Cw =G | (5c) 
3 3 
where 
A = Mass moments of inertia about principal 
axes perpendicular to the z-axis. 
Cc = Mass moment of inertia about the prin- 
Cipal axis through the hole (z-axis). 
W ,W ,W = Angular accelerations about the body 
1 # 3 fixed axes x, y, z respectively. 
WwW ,W ,W = Angular velocities about the body fixed 
1 2 3 
axes. 
GG 1G, = Moments acting about the body fixed axes. 


Since a body of revolution has only two distinct prin- 
Cipal axes, the mass moments of inertia about the x and y 
axeS are equal. The moments of inertia of the ball were 
determined by first finding the inertia tensor. From 
Refs. 21 and 22 the inertia tensor is found by starting with 
the formulation for a hollow cylinder of inside radius r and 
length 2h (Fig. 6) in which the inertia tensor in terms of 
dyadics is 


I = (U + ee)mr*/2 + (U-ee)mh*/3 (6) 


In Eq. (6) m is the mass of the cylinder (ball), U is the 
ieemdyadic and ee is the referenced dyad. Using Fig. 6 as a 


guide, the incremental value of the inertia tensor with respect 


to the center of mass, dl, is given by: 


38 






( h = /R2-9 2 ) 





Figure 6. Dimensional relationships for calculating the 
inertia tensor and moments of inertia. 
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Pie ce) p72 dm + (U - ce)h*/3 am (7) 


where 9 is the radial distance to the incremental mass, dm. 


The incremental mass is given by 
5 
dm = 41p(R*-p*) ydo (8) 


From Eq. (7), since the origin and center of mass coincide, 


I, - | ar, (9a) 


therefore, 


R R 3 
I, aI (U+ee) 20 *thydp | CUES - 4mpydo (9b) 
ig 1g 
which, after Simplification leaves 


R 3 
oe = anviveee) | 3) I fe) lays 


sa 
(9c) 

4 R S72 
$y (U-ge)/ (ae pdp 

gs 

Solving Eq. (9c) gives 
MZ 2 Be ne he 
I, = my{2(Utee) ler*+peR* (R°-r") ] + 
By 2 

$(U-ee) [Z(R¢-r?) J} - 
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Fq. (9d) is the expression for the second order inertia 
tensor of the ball obturator with outside radius R, hole 
radius r and made of a material of density y. Diadics have 
been used to describe this tensor. In order to obtain the 
moments of inertia about any axis, the diadic describing the 
tensor need only be premultiplied and postmultiplied by scalar 
multiplication by the unit vector along the axis about which 
the moment of inertia is described. It may be seen in Fig. 6 
that the e-direction vector is along the z-axis. The 
e-direction may be chosen arbitrarily, however this choice 
greatly simplifies the formulation of the integral. Thus 
the moment of inertia about the z-axis, C, 1S given by 


es | ; LL) ay ah 
Meememo sl (iitjjree)tee)] © e {[2my(er°+tzeR*(R°-r°) | 


| 1 a 
+ x(iitjjtee-ee) [mys(R*-r7) J} (10a) 


cos 


Smell tying and using the subscript s to denote particular 


values for the ball obturator: 
372 


dL 2 
c= any (=rs° 5 Ts5Fs*) (Rosas) (10b) 


Por the axes perpenducular to e through the mass center, the 


Moment of Inertia iS given by 


Me B= S-lj = atigti (11a) 
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Performing steps as in Eqs. (10) gives 


; ee ; 1 2 3/2 
A= jemi neo ree) 27 [2my (er *+yeR* (R?-r") ] 
4 B72 
ee tl Ceneey «7 [ty-=(R°-r*) ] (11b) 
Atter premultiplying and postmultiplying, 
1 2 2 2 2 a2 4 2 2 oe 
i Ty {2 [Srs + T5Rs (Rs “Is ) ] + Ty She “Is ) eave) 


In order to simplify Eqs. (5), define a non-dimensional 


X such that 


_ C-A 
A a i) 
Thus Eqs. (5) becomes 
® +Aww =G//A (5a) 
1 2 3 1 
WwW +rAww =G/A (Sb) 
2 Le 2 
w = G/C (5c) 
3 3 


The position of the ball at any time may be described by 
three angles, 8, » and wy. These are known as the Euler angles 
and their relationships to the coordinate axes already de- 
scribed may be seen in Fig. 7. Reference 20 defines Wor Ww 


2 
and w in terms of the Euler angles as 
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a = bsinOsind + bcosd (13a) 
= bsinOcosd - 6sind (13b) 
w = $ + pcosé (13¢) 


By differentiating Eqs. (13) with respect to time, w , w 
1 


and w may be obtained as 


aw = Ysin@sind + )(6cosésind + ésin@cos¢) + bcos¢ (14a) 
-8osing 

w = tsin®coss + W(Scosécosd - ¢sin@sind) - Ssind (14b) 
= 8docos¢ 

w = 6 + Ycosé - wesing (14c) 


Meesubscituting Eqs. (13) and Eqs. (14) into Eqs. (5), 
the following expressions for G Ar GSA and G 7c are found 


to be 


~sinOsingd + v(Scosésind + ésin@cosd) + §cosd - (15a) 


8osind - A(VWsinécos¢ - 6sind) (vcosé > b) = GA 


YsinOcoso + b(8cosBcosé - dsinOsino) - Bsind - (15b) 


6dcosd - A(wsinOsind + 6cosd) (Wcosé ~ b) = G/A 
Ycosé - bésine + 6 = G/e (15c) 


Transforming this system of differential equations into the 


Matrix form [A] {x} = [B] leaves 
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Sin8sind 0 cos} |v G /A - v (8cosO@sind + dsin@cos$) 


sin@cosd 0O -sind!< 6 


GA - wv (8cos8cosdé = dsin@sind) 


cos@ if 0 | 5 6 /c + wWAsine 


+ 8osind - A(Wsinécosd - 6sing) (ycosé a b) 


+ @ocosd + A(vsinésind + bcos) (Wcosé6 + 9) (16) 


By premultiplying both sides of Eq. (16) by the inverse 


of the coefficient matrix, we are left with 





. sin cos ae 

p sale —- 0 GA - p(8cos6sinod + dsin8Scosd) 
_ | Singdcosé - 5 oo . 

d¢\= Hse COs -OS2 1 G/A - (8cos&coso - dsiné@sin®d) 

5 cos -sind 0 G/e * wesind 


+ Oosinod - A(psinécosd - 6sind) (bcosé + $) 


+ Bdcosd + A (Wsinésind + 8cos¢) (Wcoseé + b) la) 


J 


In order to solve this system of non-linear second order, 
differential equations for 8, $¢ and wt, the following must be 
known: 

Initial Values (denoted by subscript 0) of 6, 6, jy, 6, b 


and w and the quantities, A, C, Ge a and x 
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Due to symmetry, 0° < 6 < 90°, and since we are interested 
in the performance of the ball away from Or? oe > Ope solnce 


the ball is attached to the projectile initially (no slippage), 
6, = bo = 0. Vo is equal to the spin rate of the projectile 
as it exits the barrel. A and C are physical constants de- 
pendent upon geometry and material properties [Eqs. (10) and 
(11)] and therefore } is a constant. $,, and to are arbitrary 
due to symmetry and the values of Gs a and G, must be de- 
termined from knowledge of forces acting on the ball. 

The solution to the problem of the motion of an axi- 
symmetric rigid body about a fixed point with no moments 
applied is well known [20,21,22,23]. Under these conditions 


there 1S no nutation and the precession rate w~ is constant. 


Thus 


and 


ac} 
(A-C) cos8® Oe 


~<e 
| 


The computer solution method was tested by substituting the 
projectile spin velocity for w Pim ilol ana Calculating 
< O 

Mueeeenie Specified value of 0 < 6, < 90 . By inserting 


these values of 6 ) and @ into the computer model, and 


oO? 
solving for @ as a function of time, 6 was indeed found to 
Be zero. In addition, all other criteria for the test case, 


1.e., constant angular momentum and constant w and $¢ were met. 
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The problem remaining is one of modeling the forces acting 


on the ball in order to determine oe cr and G.: 


B. MODELING OF THE APPLIED TORQUES 
1. Moments Due To Viscous Effects 
Without the influence of applied forces or moments, 
the angular momentum remains constant or, from first prin- 
cipals: 


—-> 


s 

M=h = 0 ene) 

In order for the ball to nutate to the position in 
Which the holes are aligned, a moment must be applied. This 
resultant torque is due to both viscous and sliding frictional 
forces. 

The first attempt was to model the torque on the ball 
as that due to laminar Couette flow of a viscous fluid between 


two flat plates where the shear stress 1S given by 


iC. = fy (20) 


The solution of the Navier-Stokes equations for this 


flow is well known and given by 


(21) 





Where dP/dx is the pressure gradient in the direction of 


Flow and, for our calculations, is assumed to be zero. The 
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governing partial differential equation is then 





2 
q“ui = 9 (22) 


By considering the relative velocity, one plate is held 
fixed while the other has velocity of magnitude Wp-F- Wr 
is the relative angular velocity about the instantaneous 
axis of rotation and Yr is the perpendicular distance from 
the instantaneous axis to the point at which the velocity 


is desired. 


The boundary conditions are therefore given as 


(23) 


where y=0 is the surface of the ball and h is the gap width 
between ball and projectile. Integrating Eq. (22) twice and 


applying the boundary conditions Eq. (23) gives 


uly) = -w rr (l-y/h) (24) 


Thus from Eqs. (20) and (24), the viscous shear 


stress at the surface of the ball, or at y=0 becomes 


E 
aL 


Ne pe (25) 
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The force due to viscous shear on the ball is given 
by the integral over the surface area of the ball of the 


shear stress Eq. (25), or 


F. = {tot (26) 


The differential area dA is given by 
dA = 2TR, (R,Sinada) (27) 


which is shown in Fig. ge. 
Siescteutina 2s.) (25) and (27) into Eq. (26) gives 
Za R. - 
TU W 
Fy, | eo ae | rsinada (28) 


h 
O 


The value of © changes from 0 at the axis of rotation 


Zoek, at a=90°. Thus 


r= R Sina (29) 


Substituting this expression into Eq. (28) 


T 
27u.R ° 
F = Hers Or J sinteas (30) 
MI 7 as 
O 


The torque imparted to the ball is the force given by 

Eq. (30) multiplied by the moment arm at which it acts. 
Again, the moment arm is a function of R, and a and is given 
Meetesing. Substituting this into Eq. (30) gives an ex- 


pression for the torque acting on the ball obturator due to 


AQ 





dA 
ts 
a 
1g 
Figure 8. Differential area, dA, orientation with respect 


to the relative angular velocity vector, w, 
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the velocity of the ball relative to the projectile ina 
viscous medium, or 


emp un 
Q= f+ [sintacs (3) 
O 
Carrying out this integration gives 
& 
TR, WHE 


- 8 _s ret 32 
ile : Siz ) 


For a given geometry and constant fluid properties, 


W) is the only variable in Eq. (32), therefore that portion 


of the expression which is multiplied by w, will be called 





ie 
the shear constant, S, such that 
= = 
Q= Sw, (33) 
where 
R 
= pe (34) 
> | lew wee: 
From the rigid body motion analysis, the relative 
angular velocity is easily determined. It is 
5 
= 2 2 2 
Wy = —W + - + - 
r Lie py) oo “ape ie aoe 2) 
where w. , W and w are the components of the projectile 


Pi Po BP 3 
angular velocity about the axes fixed in the ball. These 


relationships are given in Thomson [20] 


51 





a, = w Singsing (36a) 
Wo = w,cososing (36b) 
W = w_cos8 3 

p3 p oc) 


Equations (36) are simplified by the fact that Wr 
has only a component in the Z-direction of the inertial 
frame of reference. 

The moments Gu G. and G, are obtained by applying 
Eq. (33) and resolving the moment vector into components 
about the x, y and z-axes. The motion of the ball may then 
be described fully by Eq. (17). 


Figure 9 shows solutions for 8 as a function of time, 


using the Couette flow model with the following inputs: 


R. mee oe mre (O. 5) in. } 

fe 20 = «4.7625 mm (0.1875 in.) 

h = 0.0254 mm (0.001 in.) 

le) COOK) Ti (Ol 262 lbm/in’) (steel) 
i, = 1.92 10°° Pa-sec (4x107’71lbf-sec/ft?) 
ae = 12042.77 rad/sec (1,15x10° RPM) 


It appears that viscous effects alone are not suffi-~ 
clent to cause the ball to nutate toward @=0 in a rapid 
fashion. To check the model again, progressively larger 
Orders of magnitude of S were input into the model. The 


effect, as expected, was to cause 8 to approach zero more 
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Figure 10. 


time for applied 


6 vs. 
torques due to fluid friction 


Figure 9. 
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of Couette model (S 





BeapL~aly (see Fig. 10). From firings conducted at NWC, China 
Lake, on an indoor range, it 1s known that 8 reduces to a 
value which allows light from behind the projectile to be 
seen through the projectile by a camera aiming at the pro- 
jectile image in a mirror down range. The length of the 
range is such that, for this to occur, 8 must reach the value 
B, prior to impact on target, or in about 0.1 seconds. (See 
meee i) for Be). The viscous model described above does not 
produce the moments required to cause sufficient reduction of 
Sewithin this time of flight. 

The viscous torques on the obturator were calculated 
above assuming a solid sphere. The actual obturator has an 
outside surface area less than that of a solid sphere of the 
same radius. The surface under the influence of viscous 
shearing stress is a very complicated function dependent upon 
the angle @ as shown by the cross-hatched portion of Fig. ll. 
Calculations of the resulting moment are further complicated 
by the relative effect of the missing areas (due to the 
presence of holes) and their positions relative to the 
instantaneous axis of rotation. The time-dependent angles 
between the instantaneous axis of rotation and the z and 
Z-axes must be known in order to determine the limits of 
integration for correct moment arm consideration. 

Because the moment due to viscous effects is dependent 
upon the surface area involved and because. the surface area 


in the first calculation (assuming a solid sphere) is larger, 
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it is to be expected that by taking the holes into account 
the true moment due to viscous effects will be reduced. The 
effects of these area variations will be "blurred" by the 
high-speed relative motion and the additional complications 
necessary to account for them are not warranted in light of 
the goals of this analysis. Therefore for further develop- 
ment of viscous effects, the surface area of the obturator 
will be considered that of a solid sphere. 

Nakabayashi [24] has conducted experiments to de- 
termine the viscous torque on rotating concentric spheres 
with various fluids and gap widths. He gives the general 
expression for the frictional moment on a sphere due to flow 


in the gap as 


M = GrpReo (37) 


where 9 is the density of the fluid in the gap. For laminar 


flow 





8 -1 
ea on z Re (38) 


and for the turbulent flow 
Rome Oie2 5 


eer), O Some (774) c 


where € = h/R,. 


The Reynolds number is defined as 


2 
i 
S B 


Rens (40) 





v 





in which v is the kinematic viscosity of the fluid in the 
Gap. For the narrow gap width assumed for the BOT, the 


transition Reynolds number may be approximated by 


Re, = 70671-6 (41) 


c 


(from Fig. 2a of Ref. 24.) In the case of the BOT, the 
relative angular velocities are insufficient to assume the 
existence of turbulent flow in the gap. Thus the frictional 
moment due to viscous flow in the gap is given by 


de sn 2 eR, Be son (42) 


For small values of e€ (as is the case here), Eqs. (42) and 
(33) are identical. 
me Moments Due To Sliding Friction 

It has been noted that the extreme propellant gas 
pressures inside the barrel are assumed sufficient to hold 
the obturator fixed to the projectile during launch. How- 
ever upon exiting the barrel, the pressure distribution on 
the projectile is due to the aerodynamic characteristics of 
the projectile. Assuming the obturator is still blocking 
the projectile duct, a detached bow shock is set up (see 
Fig. 2b). The pressure, Pu downstream of this shock is given 
by one-dimensional shock theory from Eq. (1). It is assumed 
that P is sufficiently large to cause the obturator to con- 


tact the rear of the projectile socket housing the obturator. 


a 





The forces involved depend upon the projected area 
of the hole through the projectile and the contact point(s) 
of the obturator and projectile. The force acting on the 


obturator 1s assumed to be given by 
F = nr? AP [see Eqs. (1) and (3) ] (43) 


where rp 1s the radius of the hole through the projectile 
and 


AP = P -=- Pp (44) 


From sliding friction theory, the contact area is 
not a consideration; only the magnitude of the force normal 
to the contact surface and the coefficient of friction of 
the materials in contact are important. In the case of the 
obturator, the pressure force, Fy acts along the Z-axis. 
The reaction force may be broken into a component tangent 
to the obturator and a normal component acting through the 
center of the obturator. The normal component, Po! multi- 
plied by the coefficient of friction Us» gives the sliding 


Peeetion force acting to oppose the motion of the obturator: 


a =e tl (45a) 
where 

E. = F,cos8 . (45b) 
ee that 

eS FCosB Ug (46) 
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This assumes the holes through the obturator and projectile 
have the same diameter and that the point or points of contact 
are on the locus of points described by the intersection of 
the hole through the projectile and the spherical socket 
housing the obturator in the projectile. This locus of points 
is a circle of radius Ip and each point is a distance R, from 
the center of the obturator. Gravity forces are neglected 
Since they are relatively small and act both on the obturator 
and the projectile. 

For a “perfect” fit between ball and projectile, 
Pyeesliding friction force is distributed over an infinite 
number of points and, in the limit, the force at each point 
vanishes. This situation has no significance in reality and, 
therefore, engineering judgement dictates making some assump- 
tions with regard to the point(s) of application of the 
merpetion force ae 

For the general case, a point of contact P may be 
described relative to the inertial reference frame as being 
at an angle 8, from the inertial Z-axis, a distance R, from 
the center of the obturator and at an angle y from the 
X-Z plane measured toward the positive Y-axis as shown in 
Fig. 12. For this approximation, assume two contact points 


in the X-Z plane, a and P,- iB and EG are 180° ababiwon 


Biemectrcle of points so that y=0 and B=B.. 
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From the vector analysis, the projection of one vector 


onto another is given by the scalar or dot product. The 
> > 
vectors PY and P are given by: 


P ae (singi + cosgk) 
a (47) 
P = Ro Ceeisbigl{sich as cosBk) 


Their projections onto the relative angular velocity 


-> 
vector w,, with components w,, WwW, and Wz, are 


y 
> => 
P Wy = R. (w,sinB oa wcosB) 
(48) 
> > 
Pr ~ See (-w sing + w cos) 
The scalar dot product is also defined as 
-- a 
P Wy = P wycosn (49) 


> = 
where n is the angle between the vectors a and wy;. Therefore 


from Eqs. (48) and (49) 


R (w sing + w cos8) = P w cosn 
Sx Z } 5 
(50) 


R, (-w,sin8 ~ wW,cosB) = P wycose 


= > 
The magnitude of the vectors P and P is R,, so solving 


for n and 2 gives 
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nN = cos — sin8 + —  cos8 
r r 
51) 
Ww Ww 
2 = cos? Shee sing + 2. cos 
r W. 


The moment arms x and r, from Fig. 12 are therefore given by 


Ww WwW 
R,sin [cos”* (“sing + —cos8) ] 


si = 
t r r 
W W, G52) 
r =R, {1 - (sing + —cosf)?] 
1 S U) Ww 
r r 
: x vw, 
a= R.sin [cos (- psiné + pcosB) } 
(Sas), 
W 
r =R, [1 - (2cos8 - sing) ?] 
2 W- Ww 


The moment acting on the obturator due to the sliding 


friction due to aerodynamic forces are given by 


=> 
= Po Y) 
M = (xr + vane etna 
aa ( x) 5 cosBu. a (54) 


Se@estituting from Eqs. (52) and (53) into Eq. (54) leaves 


+ 


+ W W) Ww 
M = Y2p F, cospu {l - (=) 2sin?@ - (2) 2cos?28]7 cS 
ct 2 85 S W,. W, ti 
(55) 
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It should be noted that if there were n contact points that 
the normal force at each point would be given by F,/n. 
Equation (54) indicates the case of two points of contact 
discussed here. | 

As yet, there has been no attempt in the laboratory 
to produce the aerodynamic forces present on the obturator 
Beeean actual projectile in flight. The sliding friction 
forces present on the obturator in the experiments (deScribed 
in section III) are due to gravity alone. The obturator is 
Supported by the projectile simulator (See Fig. 24). Because 
the obturator is smaller than the housing in the projectile, 
a point contact at the bottom of the housing is assumed. 
Meese pOint has inertial coordinates of (-R,, 0, 0). The 


vector from the origin to the point is given by 


Bae —Roi (56) 
Using a development similar to that for obtaining the moment 
arm for the friction forces due to aerodynamic considerations, 
the forces on the model due to gravity may be determined. 
Again, the dot product is given by 

> + 


Po: Wr = -RoWy (57) 


which is the same as Pw cosé ‘ 


Therefore the angle (&) between the angular velocity 


+ < 
vector w, and Po is obtained from 


6G 





Ww 
pe="cos-. (- —* (58) 


The moment arm rg is found as in Eq. (55) from 


Ky 
ll 


R.sin [cos ? (- a) (59) 
or 


7 Wy 5 45 
25) = Re tal oe a) ] (60) 


which means the sliding friction moment due to the weight W is 


> 

> u) W 
Meee ewe. [b= (25) 2] — (61) 
Sg ss W Wr 

The torques acting on the obturator are the 
Summation of the torques due to friction, both viscous and 
Sliding (aerodynamic or gravitational). The general form 
of the torques acting to oppose the motion of the obturator 


> > 
in terms of total moment M, and the unit vector 1. along W is 
> 
G = -M,1 (62) 


Which, when broken into the various contributors becomes, 


for the projectile in flight, 


~ > 
Gp See eM.) 1, (63) 
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AERODYNAMIC MODEL 


TABLE I 


COMMON VALUES 


87.25° 


0 
0 

7669.6 Kg/m° 
7.899 mm 

4.72 mm 
3.125x107’ N-m-s? 
2.216x10~”’ N-m-s? 
0.41 

1.917x10"° Pa-s 
oa35 

SBeOSx10.° m 


12.38x10 °? N-m-s 


9.96 mm W 
76.2 mm E(8,,7/2) 
0.141 Kg re 


8.58x10 ° N-m-s? 
3.0 

1.4 

0.099 MPa 


12042 sec } 
65 


GRAVITATIONAL MODEL 


O79 37 Nh 


1.0063 


200 sec ! 





and for the experimental model, 


> > 
GC = ~ (MFM, ,) Ly (64) 


These may easily be transformed into the body-fixed coor- 
dinate system for solution of the equations of motion. 

For solving Eq. (17), a numerical integration sub- 
routine [25] was used and computer plots of @ vs. time were 
obtained. A copy of the calling program is included in 
Appendix A. Figures 13 and 14 show plotted output for the 


input parameters shown in Table l. 


C. DEVELOPMENT OF AN APPROXIMATE SOLUTION 

- Observations of the nature of the exact solution (Figs. 
13 and 14) indicate that the response of the obturator orien- 
tation (68) to the applied torques (determined from Table l 
parameters) is Similar to that of a linear damped second-order 
system with a combined step and ramp input. That is, the 


relationship closely follows a relationship of the form 


S549) (SIGNS ca = ee % 
@ = 6>tD+ Ry t-27) + DE 2 2 {[2sin[w (1-7?) ¢ + 6_] 
oO Ta (1-72) % Wn n R 


- D Siac tt-c%) “t + dnl} (65) 


in which the parameters are: 


D = step height (negative) 
R = ramp slope (negative) 
W, = natural frequency 

C = damping ratio 


dp, dp = phase angles 
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These may be determined by means of a solution to the 
linearized equations of motion under the assumption of 
small applied torques [(M/Cw,*) <<1]. This has been accom- 
plished* by a perturbation of the steady state solution given 
by Eq. (18) for the case of no applied torques. 

The first step in the approximation is to rewrite Eq. (18) 


to read 


Ap cose. + (X41) 6, = 0 (66) 


where Wor @ and 80 are constants. The dependent variables 


fo) 
of the general form of the governing equations (Eq. 17) are 
written as sums of the unforced values and the perturbation 


values. Thus 


(67) 


The general form for the applied torques on the obturator 
in the experimental model (gravity-driven sliding friction) 
1s given by Eq. (64) in which the moment term is composed of 
Geeeurd (M..) and a sliding friction (Msg) term. When both 


are expanded from Eqs. (42) and (61), Eq. (62) becomes 


a ws 5 ao 
fe - {Sw + poWR. a le a (68) 


2 , 
The theoretical development of the linear approximation 
is not complete. The analysis presented here therefore 
represents a preliminary report. 
69 





mre relative velocity w from Eqs. (¥3, 35 and 36) has com- 


ponents 
a (y-w,)sindsing + 8cosd 
Ww) (W-w..) sinOcos¢d - 4sind 
2x p (69) 
a = (p-w,,) cosé + 


To the first approximation, ,=w (Yow, for the case of 


P 
no applied torques). With this approximation, and the 


application of the perturbed variables, Eq. (67), 


Wir = Scos¢$. 

Ww) = -@sing 

2g Oo 
= ycoose. + 6. + 4 70 
fee FOOSSG Gt 8 0) 


The only term of zeroeth order in Eqs. (70) is bo SO 


mat, to the first approximation, 
wW Fw = 1 oo (ia) 


The modulus is taken in order to preserve the positive sense 


of w,- Under the same approximations as above, 


Wy. = ¢,sin? sind, (72) 


Saostituting Eqs. (71 and 72) into Eq. (60) for the 


Sliding friction moment arm, Irg/Rg becomes 
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eH 


For the 
me this 


average 


ala 


where T 


mation, 


and because of symmetry, the average is valid over one-quarter 


period, 


By letting u = 


r 
Rs 


3 


- [1-(sin6@,siny,) *] (73) 


purpose of the linear approximation, the mean value 


moment arm is required. Integrating to find the 


moment arm gives 


ae 


fh 
if (x /R,) dt (74) 


O 


is the period for one revolution. From the approxi- 


however, 
ies oe (75) 
Me Us 





therefore 
1T/4 
dg 
4 = ® 2 ® 2 2mt 76 
= + / (l-sin @,sin = ike dt (76) 
O 
2mt/T and du = (27/T)dt, Eq. (76) becomes 
T/2 i, 
za (1-sin’@ sin*u) du (77) 


galt 





Letting E(@.) denote a Complete Elliptic Integral of the 
Second Kind with modular angle aa Eq. (77) may now be 


written as 


iB 2 

R, > 7 BCS) pan 

Solutions to Eq. (78) are tabulated for various values of 

85° 
The expression for the applied torque vector (Eq. 68) 


may now be written as 











> 
u) 
G = -[S|¢] + u,WR, 2 E(6,)_= (79) 
| o| 
With these approximations, the applied torque terms 
appearing in Eq. (17) are given by 
G G 
a» , 2 -M ; 
z Sing + z-cos¢ = fy leat aa 
G G : 
xtcos} - =tsing = Fe = 
Mr al ool 
O 
G : : . > 
= = = [(p-w) cosé +o] =—=M ee (80) 
where M is given by 
Ee , Z 
M = S{o.| + USWR, = E(8)) (81) 


and Wo is the projectile angular velocity. 


eZ 





Equations (66), (67) and (80) are combined with Eq. (17). 
After simplification, and retaining only the first order 


(linear) terms, the results are: 








y Zz -(p,cote,) (82) 
a . : M , 
i (y,csc9,) 8 - = — ee! 
lo 
§ = a - Yosindy [$(1+A)-8AY sine, + prcose,}] (84) 
@ 


These equations are then integrated under the following 


initial conditions 


w (0) = Yo = wi, v = 0 

(0) = 0 $(0) = -$ = retlcos8 

6(0) = os 8 (0) = 0 

6(0) = 0 8(0) = 0 
Mee (82) yields 

y= —(M,cote.) 8 (35) 

Prom Eq. (83) 

b = -$, + (pcos) 6 = M, %o_ 

‘ NO | (86) 


By substituting the expressions obtained in Eqs. (85) and (86) 


into Eq. (84), the differential equation for @ becomes 


iS 





-Mé 





72 - posindgl¥,8 - (1+) (d, + ut) ] 


Al oo| 
ao, in a more familiar form, 


S) e 
~~ + 20g 4 @=D+ Rt 


Wn i a 
Pmere vw =) = Ws b>/|%| = -1 and 

—— 

2Alo5|V 
pee torA); sind, = - Asin26, 

Yo 
r= - SP) Seine, = - 2 
5 O5C 


Expansion of the expressions above, uSing Eq. 


that the damping coefficient 1s given by: 





1 ° 2 
6 = ———— [S]o_| + u WR.=E(6@,)] 
Sows 2° OT 
Or 
¢ = — [Ss + ugRw 2 41) 5 (,)] 
2AaV,G APocos®, 


and the ramp Slope (rate of decay) becomes 


74 


(81) shows 


(87) 


(88) 





Meo 


Mie integral of Eq. (87) with the initial conditions 


specified {6@(0) = 8 (0) = 0] is Eq. (65) where 
¢ = 2, = Qtan ? (seas (91) 


For large values of Cw it, the ramp decay dominates the 


response and the solution may be further approximated by 


= : RL 2 
@=0,+9= 06+ D+ — (wi t-22) 


n 


The time to nutate through a known angle 8 may be found 


ErOom 


t= —Q_—+ = 28 - _O _- 2 (92) 


where t, R, D and ®, are prescribed in Eqs. (88)-(90). 
Figure 15 is a plot of 6 vs. time comparing the results 

Of the exact solution from Fig. 14 with those of the 

approximation developed here. The physical constants and 


initial conditions are identical for the two cases. 


eeetHe EFFECT OF GEOMETRY ON OBTURATOR RESPONSE 
In the discussion to this point a primary consideration 


has been to model the external forces acting on the obturator, 


i 
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The effect of geometry has not been discussed, but is a most 
Significant contributor to the response of the obturator. 

C and A are functions of mass density and physical di- 
mensions. By varying the material composition (density) 
the response will vary. Utilizing the physical dimensions 
of the obturator in Table 1, the obturator volume may be 
calculated by subtracting the volumes of a right circular 
cylinder of radius oe length 2R,cos8. and the two spherical 
Caps at the ends of the hole from the volume of a sphere [26]. 


Thus the volume of the obturator is given by 


4 


2 
V = 3TR, °-¥T(R.-R,cosB) * (2R, ieecOchy = 2TL *R,cosB 


P 
Collecting terms gives 


cos*8 


= 2 
7 ty! 


a 2 - 
7 = 2mR,cos8 [R, (1 p 


By knowing the mass density, the weight of the obturator 
[for use in the torque expression, Eq, (61)] may be obtained. 
Therefore as the density varies, the values of C, A and W 
vary. Fig. 16 is a family of curves of @ vs. t for varying 
density with all other physical properties being those of 
Heble l. 

There is however a limiting case for the size of the 


hole. If the hole becomes too large (assuming same size 
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projectile hole), there may not be sufficient obturator 
Surface area to block the passage. The limiting relation- 


ship is 


This limit is purely theoretical and is not a feasible design 
parameter. Engineering judgement in the selection of realis- 
tic values iS required so that the projectile will meet the 
ballistics criteria and retain its integrity while undergoing 
the forces associated with the entire firing evolution from 
chambering to target impact. 

The effect on the obturator response due to increasing the 
hole diameter is shown in Fig. 17. The inertial contribution 
to the response (coupled term) 1s increased by increasing C/A 


and therefore A; however the friction torque (M COREL DUETOn 


sq) 


is lessened by the reduction in the mass of the ball. 
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ITI. EXPERIMENTATION 


A. EXPERIMENTAL APPARATUS 

An experimental system was designed to simulate the 
Spinning BOT. The apparatus consisted of a compressed-air- 
driven spin-up rig, optical timing mechanism, air manifold 
and associated piping, tubing and electronics. 

The obturator was fabricated from a Standard 5/8-in. 
diameter chrome~-steel (52100) bearing ball. The ball was 
annealed to allow machining and then bored along a diametrical 
axis (Fig. 18). The bored ball was then mounted in a bakelite 
metalographic specimen mount. The mounted ball was placed ina 
milling machine and a flat was machined in the bakelite at a 
Specified angle relative to the axis through the hole in the 
ball (Figs. 19 and 20). This flat was then used as the 
polishing plane for metalographic specimen preparation. A 
small flat spot was polished on the obturator to provide a 
highly reflective surface at a known orientation relative to 
to the z-axis of the obturator (Fig. 21). 

After removal from the bakelite mount, the obturator 
was placed in a three-piece, lucite housing (Fig. 22). The 
mating ends of the two hollow inner cylinders were each 
machined with a 5/8-inch end mill to a depth of approximately 
5/16-inch. When mated, a spherical cavity was formed to 


accomodate the obturator. The third cylinder was press fit 
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Figure 22. Obturator and components of lucite obturator 


housing with mating ends of inner cylinders 
machined to house the obturator. 





Figure 23. Assembled obturator housing with aluminum 
end pleces. 
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over the others after the obturator was inserted to insure 
alignment of the inner pieces and provide rigidity (Fig. 23). 
The ends of this lucite assembly were then press fit into 
aluminum end-pieces similar to those in Fig. 23. These end- 
Pieces served as the shaft for the bearings and one also served 
as the prime mover (bucket wheel) for the apparatus. The shaft 
rotated in two ball bearings mounted in aluminum pillow blocks 
aligned on a rigid pedestal. The prime mover was a bucket 
wheel machined from a solid aluminim disk (Fig. 24). 

The bucket wheel was driven by compressed air supplied 
from an installed system through an air filter to a manifold 
and then through one of two Model-10 Kendall pressure regu- 
lators. From the regulator, the air passed through a flexible 
tube to a tee, each leg of which supplied a nozzle. These 
nozzles were mounted opposite one another on the pillow block 
in such a way as to allow the air jet to impinge upon the bucket 
Wheel to cause rotation (Fig. 25). The speed of rotation was 
sensed by a Bentley Nevada Proximitor, Model 3100N, which was 
mounted above the bucket wheel to detect the passage of each 
point on the wheel. The sensor was supplied from a LAMBDA 
Regulator Power Supply Model LP413FM by 18VDC. The pulses 
generated by the sensor were counted, averaged and displayed 
as a frequency by a Monsanto Programmable Counter-Timer Model 
OB . 

While the spin-up rig was being brought up to the desired 


Speed, an air jet from a nozzle mounted rigidly on the pedestal 
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at the end opposite the bucket wheel held the obturator 

fixed to the spinning rig. This jet was supplied from the 
common air manifold and the other pressure regulator also 
through flexible tubing. The air jet passed through a hole 

in the plunger of a Rocker Solenoid, R.S. No. 10-207. This 
hole was aligned with the hole through the "projectile" to 
allow the jet to impinge upon the obturator (Fig. 26). When 
the switch at the far right of Fig. 26 was activated, standard 
L15VAC was applied to the Rocker solenoid and removed from a 
normally closed ASCO Solenoid Valve in the air supply line to 
the jet. The plunger retracted from the position shown in 

Fig. 26 and the solenoid valve closed. The retraction of the 
plunger performed three functions. It first caused a pulse 

to be generated by another Proximitor mounted next to the 
nozzle and powered from the same 18VDC supply. This pulse 
Started the timer function of a second Model 110B Counter-Timer. 
The plunger also covered the nozzle outlet by misaligning its 
hole with that of the nozzle. This removed the restraining in- 
fluence of the air jet on the obturator and prevented any air 
remaining in the supply line from impinging on the obturator. 
And lastly, in the fully retracted position, the machined and 
polished end of the plunger was positioned opposite the pro- 
jJectile hole (Fig. 27). When the released obturator nutated 
through enough of an angle, the beam of a SPECTRA PHYSICS MODEL 
132 LASER, MODEL NO. 3187, passed through the projectile and 


was reflected by the polished end of the plunger into a light 
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sensitive diode (Fig. 28). The pulse emitted by the diode 
triggered the stop channel m the timer and the elapsed time 
was displayed. Figure 29 shows the entire experimental set-up 


with air flasks in the background. 


B. EXPERIMENTAL PROCEDURE 

The obturator was positioned in the spin-up rig in such 
a way that the laser beam was reflected off the polished 
Spot back to the target mounted on the laser (Figs. 29 and 
30). The center of the target is the location of the beam 
and the obturator was adjusted to place the reflected spot 
on the center of the target. Thus the orientation of the 
z-axis through the obturator was known relative to the Z-axis 
(laser beam) through the projectile. This is the initial 
value 9,. 

The solenoid plunger was checked in position with the 
hole aligned with the nozzle (Fig. 26) and then air was 
admitted to hold the obturator in place. The pressure was 
then increased to the turbine to set the apparatus in motion. 
Figure 31 illustrates the circle traced by the reflected spot 
(Fig. 30) on the target. This circle facilitates alignment 
of the laser to insure the beam was centered through the pro- 
jectile. (It can be seen in Fig. 31 that the beam was slightly 
Off-center.) Alignment was accomplished by positioning the 
beam from side to side for horizontal alignment, then raising 


Or lowering the spin-up rig's mounting platform with adjusting 
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Figure 30. Laser reflection on the target from the 
polished spot on the obturator. 





Figure 31. Circle traced by the reflected laser while 
the apparatus is Spinning. 
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screws for vertical alignment. When aligned, the circular 
trace of the reflection exactly follows the circles of the 
target. 

Once the alignment had been checked and the obturator 
positioned at the known OO! the pressure wasS increased to 
the turbine nozzles until the desired spin rate was achieved. 
The regulator was then adjusted so as to maintain the spin 
rate at +5 Hz on the digital display or approximately +1.75 
rad/sec. (The actual tubular projectile spin rate is 
approximately 1.2x10"* rad/sec). Once the desired rate was 
achieved and noted, and the timer checked and reset if neces- 
Sary, the switch (Fig. 26) was activated. This secured the 
air to the nozzle holding the obturator and retracted the 
solenoid plunger (Fig. 27). The obturator was released and 
began to move relative to the spin-up rig. When the angle 
8 reached the value at which the laser beam could pass through 
the obturator hole (0=8,,=8), the beam struck the polished 
plunger and was reflected into the light sensitive diode 
housing (Fig. 28). The reflected beam striking the diode 
caused the diode to emit a pulse which turned off the timer 
gate. The elapsed time from plunger retraction to release 
the obturator, to the obturator nutating to 8. was dis- 
played on the timer display to the 0.0001 second. This time, 


ic was recorded. 
a= ecorded 


o) 





oe DATA REDUCTION 

The only data reduction required was the conversion from 
Hz to Radians per second. The frequency read-out on the 
digital display indicated the number of teeth on the bucket- 
wheel pasSing the sensor in one second. Since there were 
eighteen teeth on the bucket wheel, the conversion was 


simply 








Appendix B is a listing of the raw data and the 


corresponding values of Wa 
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ie COnRRGUALLON OF RESULIS 


Figures 32 through 37 show 6 vs. time for various values 
of wu, Peomet0 to 1000 sec’*’. At high and low values of Wor 
§ approaches zero more slowly than at the intermediate values. 
Figure 38 is a plot of Eq. (92) and the experimental data. 
The plot of Eq. (92) uses the physical qualities from Table l 
and u.=0.35 over the range 100<w,,£1000 Seclo) it ds clear 
that the experimental data falls below the approximation curve 
at high We and rises sharply above the curve at low We ° 

In order to determine what parameters may effect the 
position of the approximate curve, sSenSitivity checks were 
run on the quantities in which the most uncertainty existed, 
namely Ug and ®, (even though @ could only vary by #1°). 
Figures 39 and 40 indicate the sensitivity of tor to variations 
ilygl Ue and @, respectively. Small variations in oe are not 
discernable within the range of scatter of the data (Fig. 41). 
Variations in u, have a much more pronounced effect. From 
Fig. 42 it may be seen that the data at higher values of Ws 
correlate well fora U=0.95. 

It has been noted that the approximate solution is based 


upon the assumption M/Cw,* 


<<l. At values of w, (<200 S867"), 
this assumption is no longer valid and therefore calculations 
and correlation of the data to the approximate solution will 
not be conducted for Ww, 200 Sccmeetsomat low values of Wor 
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-_ appears to approach » almost assymptotically. Because 
652909 (87.259), the obturator may be near a metastable state 
at low Wy such that the coupled terms in the equations of 
motion have little effect. The amplitude of oscillation at 
low Bp is also relatively small. 

meegeneral, for large values of ©4.t (Eq. 65) and 
neglecting the effect of viscosity, t becomes 


Av) 
_ Atl M p ae 
t = T= ———__ + __#- [(6,-8) ssind)] (9) 33), 


Awp ?coS 8, Msin§, 


where M is the gravitational contribution from Eq. (81). 


This may be written as 


t = 1 +7 wy (94) 





with T and T being constants determined from the known 
physical parameters (Table 1). Utilizing values of 

Gee 87.259, 9 = 0. = sin™’ (r/R,) = 36.734° and the measured 
weight of the obturator, W = 0.0798 N, 

Mis found to be 


M = 6.34x107* LD N-m 


Therefore, 


| 
ll 


2.039x10° 


and 


rH 
ll 


3.03xlo * / 1 (95) 
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thus 


0 2.030 


eS 


u t+ 


But for large values of w_, 


Pp 


3.03x107*w 


p 
v 


B20sx10" 
D 


(96) 


A least squares fit of the data (Wp7235.97) Glo) HCl G AGA) F 


such that t = T Wp, gives 


2 Uti Woi 


Substituting in 


2 
Dts 


2 
LW 4 


2 2 
Lty ea 


therefore 


T = 8.6603508x10 


the data points (denoted by 


2.1046884 sec? 


2.8413714x10°sec ? 


2.430257x10? 


4 


2 


sec 


txt 


Since the standard deviation for T is given by 
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in Appendix B), 








where 


the following results are obtained: 


iS 
il 


fe) 1.4659547x10? 
On? = 3.9181053x10 7° 
= th, Cala) Be 

‘o} 

=e = 0)-023 

iE 


The data therefore, deviates from the least squares value 
by only 2.3 percent. 


Using the value calculated for T rec eo), 


= h 
gy = 3203x107" _ 9.35 
8.66x10 


meth £(87.25°) = 1.00458, 


Yr = 2 (1.00448) = 0.63954 
fe = 7 
S 
thus 
0.35 
Us = 0eO395T = 0.5473 


FEAISE 





With this value of Mos values of ee were obtained from 
the exact solution over the range 125 <wy</00 sec’ and 


plotted with the approximation of t.,, from Eq. 94 and the data 


a 


Meents (Fig. 43). 
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V. CONCLUSIONS 


The response of the obturator (8) as a function of time 
to applied torques was found to be that of a rigid body of 
revolution about its mass center. Euler's Modified Equations 
of Motion in terms of Euler Angles describe the response to 
these applied torques. The response is predictable and was 
compared for various combinations of physical parameters. By 
varying density, hole size and angular velocity, various re- 
sponse times may be obtained. It was found that in each 
case, a minimum response time was attainable for a given value 
of Wp For the case of the experimental apparatus, this mini- 
mum occured at a ast sec !, 

The analytical model and a linear approximation to this 
model were used with the parameters of the experimental 
apparatus and close agreement to the data was achieved at 
Wr > 200 sec '. Below Wy = 200 sec~! the exact solution appeared 
to be higher and the approximate solution lower than the experi- 
mentally observed values of toy. This is explainable in the 
approximation since it was based upon the assumption that 
M/Wp°C<<l. At low values of wp, this assumption does not hold. 

The trend in the exact solution at low w, is to approach 


ead 


tor = © almost assymptotically. This trend is observable in 
the experimental data as well and may be explained by the 


metastable initial position of the obturator @,=90° (8)=90° is 
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a stable position). The effect of the coupled terms in the 
equations of motion is small at low ays therefore the initial 
effects are not as great and the amplitude of oscillation of 
8 1s insufficient to cause departure from this metastable 
position. Since the actual projectile operates at Wpy=12,000 


sec-!, high values of Wy are of interest. The theory agrees 


1 and therefore the 


well with the data above Wp=200 sec. 
correlation to the data was made at these higher values. 

Because the linear approximation agreed closely at higher 
Wor a least squares fit of that approximate equation to the 
experimental data (Wp*235.97 Se cmeleacconthe calculatiton 
of a higher sliding friction coefficient than had been used 
previously. The assumed value was He-0 +2 for laminated 
plastic on steel [27]. This was found to be too low by the 
least squares fit. A sensitivity check of the approximate 
Bemutiton for t,, tO variations in UY, also indicated that 0.35 
was too low (Fig. 39). Therefore, a method for determining 
the sliding friction coefficient from the experimental data 
and the approximate solution was discovered. 

By knowing initial orientation (8,), projectile spin rate 
(Wy) and the physical dimensions and material properties of 
the projectile, the response may be predicted in a gravity 
environment with reasonable accuracy. 


The aerodynamic model appears to give reasonable results 


within the scope of this study. The complex forces acting on 


NS) 





the obturator in flight are extremely difficult to model 
theoretically and therefore they were simplified greatly. 
The confidence gained in the gravity model, however, leads 
to the conclusion that once the complex aerodynamic forces 
are better understood, the motion of the ball within the 


flight model will be predictable. 
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VI. RECOMMENDATIONS 


Further experimentation should be undertaken with 
varlation of other parameters (hole size, obturator material) 
to obtain more confidence in the analytical model. Wind 
tunnel tests of the BOT should be conducted in order to more 
closely approximate the actual forces on the obturator in 
flight, and to determine the position (60) at which the de- 
tached bow shock is swallowed. This would greatly enhance 
the predictive capabilities of the model and thus allow 
coupling to a trajectory model for predicting flight char- 
acteristics. This, in turn, could be coupled to a numerical 
optimization routine for optimizing the BOT design for the 


desired performance criteria. 
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APPENDIX B - RAW DATA 


READOUT TIME W READOUT TIME We 
[Hz ] [sec]. [sec *] [Hz] [sec] [sec ?] 
776 2 a 270.88 880 ea al Ua 307.2 
676 5 PLAS 235.97 920 - 3026 321. la" 
679 -2181 237.02*% 920 -2847 321.14% 
5/73 feos? )6200.01 936 -2816 326.7% 
70 P2079 198.97 983 - 2699 343.1* 
574 -2009 200.36 LOOL - 306 349.4% 
570 - 197 198.97 994 -2972 346.97* 
Sal -2031 199.32 3055 - 296 368.26% 
475 -2099 165.8 1063 nes 18) 371.06* 
479 wives s 167.2 aso Peels) 401.77* 
5/70 -2179 198.97 1022 . 369 356.7 * 
475 25 165% ys 378 410.15% 
479 2223 L6G 4:2 L253 - 3348 439.12*% 
470 ~ 230 164.1 1278 - 3968 446.1 *® 
428 -1831 149.4 1255 -4268 438.08x% 
a25 -1786 148.35 Pa al - 3088 426.2 * 
430 mtGooo 150.1 1170 -284 408.4 * 
367 ihegio 128.1 1158 52 Gee 404.2% 
363 eeauo 7? 6s L 26. 7 
370 -1959 129.15 
Ot | 2249 96.7 
267 2292 93.2 
279 2259 97.4 
323 - 230 112.7 
B23 mao” V2. 7 
Ds Po 493 To 25 
22 3 658 TIes8 
223 599 77.8 
226 S241 78.9 
807 2592 281.7% 
810 -2339 282.7% 
853 2773 297.8% 
853 -26031 297.8% 
850 -2917 296.7% 
854 ~267 298.1* 
921 - 3022 # 321.49 
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